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Abstract 

An algorithm is presented to fit precipitates in atom probe tomographic data sets 
as equivalent ellipsoids. Unlike previous techniques, which measure only the radius 
of gyration, these ellipsoids retain the moments of inertia and principle axes of 
the original precipitate, preserving crystallographic orientational information. The 
algorithm is applied to study interconnected 7'-precipitates (LI2) in the 7-matrix 
(FCC) of a Ni-Al-Cr alloy. The precipitates are found to coagulate along (llO)-type 
directions. 
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1 One of the principle challenges that analyzing three-dimensional atom-probe 

2 tomographic (APT) results poses is the amount of raw data that the instru- 

3 ments are now able to collect [1,2]; we have collected continuous data sets as 

4 large as 2.1 x 10 8 properly ranged atoms from a single specimen. It is neces- 

5 sary to extract information about spatial and compositional measures, such 
e as precipitate size from these data sets. A common method to measure the 
7 size of precipitates is to calculate a single radius of gyration (or, from this, the 
s Guinier radius) [3,4]. While this technique works well for equiaxed, spheroidal 

9 precipitates [5,6], many alloy systems investigated by APT have non-spherical 

10 features [7,8,9,10,11,12,13,14,15]. Improperly reconstructed data sets, which 

11 have preferential evaporation or local-magnification effects, may also exhibit 

12 non-spherical features [16]. Techniques that measure the angular eccentricity 

13 of features can be used to quantify the accuracy of a reconstruction. The radius 

14 of gyration technique does not, however, retain three-dimensional information 

15 concerning precipitate orientation. 

16 In this article, a more general alternative to the radius of gyration is presented 

17 and applied. Best-fit ellipsoids have equivalent centroids, moments of inertia, 
is and principle axes for arbitrarily shaped precipitates. The crystallographic 

19 orientations of the resulting ellipsoids are then used to study the coagulation- 

20 coalescence coarsening mechanism in a Ni-Al-Cr alloy, which occurs when 

21 the 7'-precipitate number density is large (> 10 24 m -3 ) and the edge-to-edge 

22 distance between adjacent 7'-precipitates is small (< 2 nm) [17,18]. 

23 In lattice kinetic Monte Carlo simulations, a coagulation-coalescence coarsen- 

24 ing mechanism is reported [18]. This mechanism is caused by non-equilibrium 

25 overlapping diffusion fields, which originate from the long-range vacancy-solute 

26 binding energies and a small mean edge-to-edge interprecipitate distance. The 

27 non-equilibrium concentration profiles observed at the 7'-precipitate/7-matrix 

28 interfaces lead to a higher interfacial free energy than for fully equilibrated 

29 7'-precipitates. The excess free energy of the region of overlapping concen- 

30 tration profiles ("diffuse neck") can decrease by changing the concentration 

31 thereof into a well-formed neck [18]. Phase-field simulations find that the rate 

32 of 7'-precipitate coalescence is increased when the 7/ / 7 / -interfacial width is 

33 increased artificially, thereby increasing the overlapping diffusion fields [19]. 

34 While nanometer-sized coagulated 7'-precipitates might have been observed 

35 experimentally, past studies only commented on whether precipitates appeared 

36 to be non-equiaxed [20,21,22], necked [23,24], or chemically ordered [17]. They 

37 did not explore the crystallographic orientation for precipitate coagulation. 

38 A Ni-5.2 Al-14.2 Cr (at.%) alloy was melted under an Ar atmosphere and 

39 chill cast. Its chemical composition was verified by inductively coupled plasma 

40 spectroscopy. The alloy was homogenized for 24 h at 1300°C, which resulted 

41 in coarse grains (0.5-2 mm diameter). After homogenization, the alloy was 

42 annealed at 900°C (7-phase field) and water quenched to ambient tempera- 
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ture. The solutionized alloy was sectioned and aged for 4 h at 600°C and then 
quenched. This treatment leads to the greatest percentage of 7'-precipitates 
that are interconnected by necks (30 ± 4%) [17]. The specimens were cut, 
ground, and then electropolished into APT tips. Three separate APT runs of 
ca. 3 x 10 6 atoms were collected on a first-generation 3-D APT [25,26] at a 
specimen temperature of 40.0 ± 0.3 K, a pulse fraction of 19%, and a pulse 
repetition rate of 1.5 kHz. The computer programs IVAS (Imago Scientific In- 
struments) and ADAM [27] were used to analyze APT data. The 7/7' interface 
is delineated using a 9 at.% Al isoconcentration surface [28] and the atoms 
contained within the 7'-surface were exported and segmented into individual 
7'-precipitates by a modified envelope algorithm [3,4,29,30]. 

The 7'-precipitates are divided into three classes: (i) single 7'-precipitates (un- 
coagulated, without a concave neck) that are not cut by the surface of the ana- 
lyzed volume; (ii) two or more coalesced 7'-precipitates that are interconnected 
by a concave neck; (iii) 7'-precipitates cut by the analysis volume boundary. 
Class (i) accounts for 42% of 7'-precipitates analyzed, class (ii) accounts for 
28% of the 7'-precipitates analyzed, and only 12% of them are formations of 
more than two 7'-precipitates (the largest of which is made up of five dis- 
tinguishable 7'-precipitates). The best-fit ellipsoid method yields quantitative 
results for all three classes; particularly class (ii), which is important for under- 
standing the coagulation-coalescence mechanism of 7'-precipitate coarsening. 
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[Fig. 1 about here.] 

A schematic that explains the fitting of an equivalent ellipsoid to atoms of 
coagulated and coalesced precipitates is presented in Fig. 1. For a reference 
space defined by Cartesian axes Xi, X 2 , X 3 (typically the analysis direction 
and the two orthogonal principle directions of the area detector), the major 
and minor axes of the best-fit ellipsoid of a precipitate containing iV atoms 
are determined directly from its eigenvalues (Ai, A 2 , A3), also referred to as 
the principle axes [31]. The principle axes are segments along the transformed 
Cartesian X[, X%, and X^ axes and are obtained from the diagonalization of 
the characteristic length matrix, L (also known as the inertia tensor). This 
diagonalization is obtained by a Jacobian transformation [32] of a symmetric 
second-rank tensor, as follows: 
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76 where the characteristic lengths, Ijk, are calculated from the positions of i 

77 atom in the reference space (xi(i), x 2 (i), and x 3 (i)), relative to a precipitate's 
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78 center of mass (xi(com), £2 (com), and, £3 (com)), averaged over N atoms, 

79 employing: 



1 N 



l kk= t;E E foW ~ Xj(com)) (2) 

1 * 

/jfc = - j-; Yj ~ ^j(com)) (x fc (z) - x fc (com))) for j 7^ fc (3) 



so The diagonalization of the L matrix follows a procedure outlined in Ref. [32]. 

si The transformation matrix used for this diagonalization yields the orientation 

82 of the ellipsoid with respect to the reference state. Defining A x > A 2 > A3, the 

83 semi-axes (Si) of the best-fit ellipsoid are given by: 



Si = J - (A,- + A fc - A<) for j + k; (4) 



84 where S3 > S2 > Si are the major semi-axis and two minor semi-axes, respec- 

85 tively. 

se [Fig. 2 about here.] 

87 Figure 2 shows a reconstruction of one APT data set where ellipsoids have been 

as fit to interconnected 7'-precipitates. As is seen, the best-fit ellipsoid retains 3- 

89 D size and crystallographic orientation information that is jettisoned by other 

90 techniques [3,4]. 

91 Single 7'-precipitates that are not cut by the analysis volume have aspect 

92 ratios of I 2 - = 1.5 ± 0.5 and I 2 - = 1.3 ± 0.2. The closeness of these values to 

02 01 

93 unity is consistent with equiaxed, uncoagulated precipitates. 

94 Coalesced 7'-precipitates that are interconnected by necks have aspect ratios 

95 of # = 2.9 ± 0.9 and #1 = 1.3 ± 0.3. The ratio for fa is about twice the 

*2 Jl J2 

96 same ratio for uncoagulated precipitates, but |£ is about the same for the two 

97 classes. This demonstrates that a majority of these consist of two equiaxed 

98 7'-precipitates that have coagulated and undergone coalescence. 

99 Precipitates that are cut by the edge of the analysis volume can serve as a check 

100 of the best-fit ellipsoid method. These have g = 2.2 ± 0.8 and f| = 1.4 ± 0.4 

101 because the majority are equiaxed, uncoagulated 7'-precipitates that are, on 

102 average, cut in two by the analysis boundary. Sampling bias by coagulated 
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103 7 /- P rec ipitates is small, as they make up a minority of precipitates and many 

104 can be isolated (if their concave necks are in the analysis volume). 

105 [Fig. 3 about here.] 

loe It is useful to relate the axes of the analysis volume with specific crystallo- 
107 graphic directions to study the orientation of 7'-precipitate coagulation. The 
los analysis direction, which was chosen near the 002 crystallographic pole, pro- 
log vides us with the [001] direction. The [010] and [100] directions can be de- 
no duced from a Field-Ion Microscope (FIM) image (Fig. 3) [33] to within 1-4°. 
in FIM micrographs also demonstrate that local magnification effects are negli- 
n2 gible in this alloy, which is consistent with the small lattice parameter misfit, 
us 5 = 0.0006 ± 0.0004 [17). 

H4 [Fig. 4 about here.] 

H5 The rotation matrix used in the best-fit ellipsoid method yields the Bunge 

lie Euler angles [34] for the major principle axis of the 7'-precipitate with re- 

H7 spect to this crystallographically-resolved reference system. From these, an 

us inverse pole figure (Fig. 4) shows the orientations of 7'-precipitates that are 

n9 the result of the coagulation-coalescence coarsening. There is a preference for 

120 coagulation along the (llO)-type directions. 30% of the coalesced precipitates 

121 are within 10° and 71% are within 15° of (llO)-type directions. In the FCC 

122 structure of the 7-matrix, (110) is the fastest diffusion path for solute clus- 

123 ters, and is therefore consistent with the model presented in Ref. [18]. Some 

124 7'-precipitates coagulated along (100), which is the next-fastest diffusion path 

125 and none coagulated along (111), which is a slower path. 

126 No crystallographic orientational preference was found for single, uncoagu- 

127 lated 7'-precipitates. This supports, once more, the equiaxed nature of 7'- 

128 precipitates and the proper reconstruction approach (with negligible local 

129 magnification) for analyzing the raw APT data. The 7'-precipitates that are 

130 cut by the analysis boundary show a preference for the [100] direction, as that 

131 direction makes up the majority of the analysis frustum's surface. 

132 The measurement of size and orientation of non-spheroidal precipitates in APT 

133 data requires more spatial information to be preserved than for the commonly 

134 used methods currently reported to date. We have demonstrated that the 

135 best-fit ellipsoid technique preserves the center of gravity, moment of inertia, 

136 and principle axes of any precipitate. The technique is applied to specific re- 

137 suits for a Ni-Al-Cr alloy with both uncoagulated equiaxed 7'-precipitates and 

138 nonequiaxed lobed precipitates that formed through a coagulation-coalescence 

139 coarsening mechanism. 
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255 1 The best-fit ellipsoid of a precipitate is determined in three 

256 steps: (i) for a reference set of axes (X 1; X 2 , and X 3 ), identify 

257 the Xi, X2, £3 coordinates of all atoms in a precipitate and 

258 its center of mass to obtain a moment of inertia tensor 

259 (L) (Eqs. 1-3); (ii) the orientation of the principle axes 

260 (X[, X' 2 , and X' 3 ) of the ellipsoid axes construct a Jacobian 

261 transformation matrix that will diagonalize L. (iii) the lengths 

262 of the major and minor semi-axes of the best-fit ellipsoid are 

263 found from the eigenvalues of the transformed matrix. (Eq. 4). 11 

264 2 A representative data set for a Ni-5.2 Al-14.2 Cr (at.%) 

265 alloy, whose thermal history is discussed in the text. The 

266 coloring scheme matches that of Fig. 1, with Ni atoms in 

267 green, Al atoms in red, and Cr atoms in blue. The atoms in 

268 the 7-matrix (FCC) are omitted for clarity. The gold colored 

269 best-fit ellipsoids indicate those 7'(L1 2 ) precipitates that are 

270 interconnected by necks. 12 

271 3 An FIM image, centered on the 002 pole, taken before APT 

272 analysis. The small square denotes the area and orientation 

273 for the data set in Fig. 1. The larger square connects the 113 

274 family of poles and the edges of this square give the [010] and 

275 [100] directions. 13 

276 4 An inverse pole figure for all interconnected 7' (LI2) 

277 precipitates in this study, presented in the standard 

278 stereographic triangle. There is a preference for coagulation 

279 along and close to (llO)-type directions, which is consistent 

280 with a diffusion controlled mechanism for coagulation. 14 
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Fig. 1. The best-fit ellipsoid of a precipitate is determined in three steps: (i) for a 
reference set of axes (X\, X2, and X3), identify the x±, X2, X3 coordinates of all 
atoms in a precipitate and its center of mass to obtain a moment of inertia tensor 
(L) (Eqs. 1-3); (ii) the orientation of the principle axes (X[, X' 2 , and X' 3 ) of the 
ellipsoid axes construct a Jacobian transformation matrix that will diagonalize L. 
(hi) the lengths of the major and minor semi-axes of the best-fit ellipsoid are found 
from the eigenvalues of the transformed matrix. (Eq. 4). 
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Fig. 2. A representative data set for a Ni-5.2 Al-14.2 Cr (at.%) alloy, whose thermal 
history is discussed in the text. The coloring scheme matches that of Fig. 1, with Ni 
atoms in green, Al atoms in red, and Cr atoms in blue. The atoms in the 7-matrix 
(FCC) are omitted for clarity. The gold colored best-fit ellipsoids indicate those 
7 ; (Ll2) precipitates that are interconnected by necks. 
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Fig. 3. An FIM image, centered on the 002 pole, taken before APT analysis. The 
small square denotes the area and orientation for the data set in Fig. 1. The larger 
square connects the 113 family of poles and the edges of this square give the [010] 
and [100] directions. 
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Fig. 4. An inverse pole figure for all interconnected 7' (LI2) precipitates in this 
study, presented in the standard stereographic triangle. There is a preference for 
coagulation along and close to (llO)-type directions, which is consistent with a 
diffusion controlled mechanism for coagulation. 



11 



